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折扣 {0-1} 背 包 问 题 粒子 群 算法 的 贪 禁 修复 策略 探究 
代 祖 华 ， 周 XX, Em ERR 


(西北 师范 大 学 计算 机 科学 与 工程 学 院 ， 兰 州 730070) 


摘 要 : 群 智 能 启发 式 算 法 求解 折扣 {0-1} 背 包 问 题 (D{0-1}KP) 时 ， 为 提升 求解 效率 和 求解 质量 ， 需 采用 某 种 修复 与 
优化 策略 将 非 正常 编码 个 体 转换 为 符合 解约 束 条 件 的 编码 个 体 。 在 引入 项 集 价值 密度 概念 基础 上 ， 以 粒子 群 算法 
(PSO) 为 例 ， 提 出 一 组 基于 项 集 的 贪 禁 修复 与 优化 方法 (Group Greedy Repair and Optimization Algorithm, GGROA), 
并 进一步 构造 PSO-GGRDKP 算法 (PSO based GGROA for solving D{0-1}KP) 以 探 完 GGROA 方法 的 可 行 性 和 性 能 。 
PSO-NGROADKP(PSO based NGROA for solving D{0-1}KP) 和 PSO-GRDKP(PSO based GROA for solving D{0-1}KP) 
是 基于 项 贪心 修复 与 优化 方法 的 粒子 群 算法 。 在 D{0-1}KP 标准 数据 集 的 实验 结果 表明 : 5 PSO-NGROADKP 和 
PSO-GRDKP 相 比 ，PSO-GGRDKP 算法 的 解 误差 率 略 高 ， 但 算法 时 间 性 能 分 别提 升 13.8%、12.9%。 

关键 词 : 折扣 {0-1} 背包 问 题 ， 启发 式 算法 ; 粒子 群 算法 ; 非 正常 编码 个 体 ; 贪心 修复 与 优化 ; D{0-1}KP 数据 集 
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Greedy repair strategy of particle swarm optimization for discounted {0-1} knapsack problem 


Dai Zuhua, Zhou Bin, Long Yujing, Wang Zongquan 
(College of Computer Science&Engineering, Northwest Normal University Gansu 730070, China) 


Abstract: Swarm intelligence heuristic algorithm is used to solve discounted (0-1) knapsack problem (D (0-1) KP) . In order 
to improve the solution efficiency and quality, a repair and optimization strategy is needed to convert abnormal coding 
individuals into coding individuals that meet the solution constraints. On the basis of introducing the concept of group value 
density, taking particle swarm optimization algorithm (PSO) as an example, a set of greedy repair and optimization methods 
based on group (GGROA) is proposed, and the PSO based GGROA for solving D{0-1}KP algorithm (PSO-GGRDKP) is 
further constructed to explore the feasibility and performance of GGROA. PSO-NGROADKP and PSO-GRDKP are PSO 
algorithms based on item greedy repair and optimization method. The experimental results on D {0-1} KP standard data set 
show that compared with PSO-NGROADKP and PSO-GRDKP, PSO-GGRDKP has slightly higher error rate, but the time 
performance of the algorithm is improved by 13.8% and 12.9% respectively. 

Key words: discount (0-1) knapsack problem; heuristic algorithm; particle swam optimization algorithm; non-normal coding 
individual; greedy repair and optimization; D (0-1) KP dataset 
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10-1) E E [RE BI( (0-15 Knapsack Problem, {0-1}KP) 是 计算 项 集 i 的 各 项 价值 密度 记 作 oon- 4.”%.", EAA 

机 科学 一 个 重要 的 NP Complete 问题 ， 也 是 一 类 经 典 组 合 优 至 多 有 一 项 被 选择 装 入 背包 。 在 不 超过 背包 载重 量 C 的 条 件 下 


[n 
化 问题 ， 在 商业 、 经 济 、 管 理 、 安 全 等 人 领域 有 着 广泛 应 用 背 ”从 给 定 项 集 选择 满足 装 入 背包 要 求 的 项 ， 使 得 装 入 背包 所 有 项 上 
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景 。{0-1}KP 自 被 提出 后 的 几 十 年 里 被 反复 地 研究 ， 衍 生出 价值 系数 之 和 达到 最 大 。 


可 漳 法 和 分 支 限 界 法 等 ， 非 精确 算法 主要 有 随机 算法 、 近 似 x; =1(0< j<3n-) RRT J RAKE, *%=0 表 示 项 j 未 装 入 背包 ， 
算法 、 生 物 算 法 等 .2005 Æ, H GUDER 首次 提出 的 折扣 {0- ”对 于 决策 向 量 Q2, D(0-1) KP. 的 整数 规划 模型 为 
1} 背 包 问 IBi(Discount {0-1} Knap-sack Problem, D (0-1! KP)U?1 e 
max f(x) = max 5 Ga; ps; + Xia Psi + Xia Paisa) 
是 {0-1} 背 包 问题 的 拓展 形式 ， 作 为 刻画 折扣 销售 、 捆 绑 销售 a UBER () 
等 商业 活动 现象 的 经 358055 j 模 型 ， AE [ni 家 设计 商业 ET T4847 3E. Fk 中 : Xi 十 Jp Xo El Xi Xi X2 E {0,1} (2) 
H RAW 买 商 h 决策 的 和 学 计算 依据 ° » S uw, + Xi Wairi TOX443Wy3) ES C (3) 
D{0-1}KP 实例 由 一 组 项 集 item_group) 组 成 ， 每 个 项 集 有 3 Jan i0 
(项 ，item 可 供 背 包装 入 选择 ， 定 义 1 给 出 了 D(O-L)KP 第 一 数学 模型 。 D{0-1}KP 的 项 集 有 4 种 选择 状态 ， 如 果 背 包容 量 和 项 的 
定义 1 D1{0-1}KP 第 一 数学 模型 站 ,给 定 n 个 项 集 和 载 价值 系数 取 值 范围 很 大 , 则 不 选择 某 个 项 集 的 约束 条 件 较 难 确 


um 


为 C 的 背包 ， 每 个 项 集 ii=0,4…n-D 由 3 个 项 (item) 组 成 ， 定 ， 这 意味 着 D{0-1}KP 的 算法 难度 要 大 于 {0-1}KP UI, 


精确 算法 、 非 精确 算法 两 大 类 算法 。 精 确 算法 有 动态 规划 、 D{0-1}KP 是 一 个 特殊 整数 规划 问题 ， 定 义 决策 变量 


以 {0- 


编号 记 作 353i 13i 2, 项 对 应 重量 系数 记 作 war 各 项 DRP 研究 成 果 为 基础 ， 学 者 们 进一步 研究 了 D(0-1)KP 的 各 
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类 算法 。 其 中 基础 动态 规划 算法 Sl(Basic Dynamic Programming, ”发 式 算法 求解 性 能 会 有 差异 。 本 文 基于 项 集 价 值 密度 印 概 念 ， 
BDP) 是 处 理 小 规模 数据 的 精确 解 算法 。2016 年 贺 角 朝 等 人 负 构造 一 组 基于 项 集 的 贪 禁 修复 和 优化 方法 , 并 探究 其 可 行 性 和 
以 “所 选择 项 的 价值 系数 之 和 使 总 重量 最 小 ”原则 构造 动态 规 性 能 , 进一步 讨论 贪 禁 修 复方 法 与 数据 实例 类 型 之 间 的 适应 性 
划 目 标 函 数 ， 提 出 一 种 新 动态 规划 算法 (New Exact algorithm for 关系 。 为 避免 启发 式 算法 差异 对 研究 内 容 的 影响 ， 以 粒子 群 算 
D{0-1}KP, NE-DKP)， 当 背包 容量 大 于 所 有 的 项 集 第 三 项 价值 累 。 ”法 为 例 来 构造 D(0-1) KP 求解 算法 。 
加 和 时 ，NE-DKP 算法 性 能 好 于 BDP 算法 。 动 态 规划 算法 是 伪 Js 
多 项 式 时 间 复杂 度 ， 一 般 不 适用 于 大 规模 D{0-1}KP 实例 的 求解 。 1 MAMATE 

群 智能 启发 式 算法 是 近年 来 求解 D{0-1}KP 的 主流 算法 ， 1.1 二 元 粒子 群 优 化 算法 
对 于 解决 大 规模 实例 有 突出 性 能 优势 。 应 用 D{0-1}KP 第 一 数 Kennedy 和 Eberhart 通过 对 乌 群 捕食 行为 的 研究 ,在 1995 
学 模型 设计 启发 式 算法 ， 决 策 向 量 098.05 的 不 同 取 值 组 年 提出 标准 粒子 群 算 法 (Particle Swarm Optimization, PSO) [9]， 
合 对 应 不 同 个 体 , 这 种 二 元 编码 法 虽然 便于 个 体 演 化 算 子 的 实 BPSO(Binary Particle Swarm Optimization) 是 一 种 应 用 于 离散 空间 
现 ， 但 由 于 可 行 解约 束 条 件 (2)(3) 的 限制 ， 编 码 空间 中 非 正常 。 搜索 的 二 元 粒子 群 优化 算法 7，Bansal 和 Deep 用 粒子 速度 作为 
编码 个 体 ( 即 个 体 编码 不 对 应 问题 可 行 解 ) 的 概率 至 少 为 1-(3) B1, {0-1}KP 中 物品 选择 为 1 或 0 的 概率 ， 提 出 了 一 种 修正 二 元 粒子 
为 此 ,需要 采用 某 种 策略 将 非 正常 编码 个 体 转换 为 符合 解约 束 群 优化 (Modified Binary Particle Swarm Optimization，MBPSO) 算 
条 件 的 编码 个 体 , 以 提升 求解 效率 和 求解 质量 . Michalewicz[ 人 1、 法 来 解决 背包 问题 (中 。MBPSO 算法 原理 描述 为 : 假设 搜索 空间 
贺 雪 朝 B4 等 学 者 先后 提出 的 贪心 修复 与 优化 策略 是 消除 非 正 D 维 ,定义 D 维 向 量 二 = Gusta? 表示 粒子 群 中 第 ; 个 粒子 
常 编码 个 体 效 果 较 好 的 方法 ,文献 [4] 在 研究 D (0-1) KP 粒子 群 的 位 置 ， 对 应 粒子 速度 为 w=Wovow…w…%vp) ， 个 体 极 值 
求解 算法 时 , 证 明了 给 定 项 集中 三 个 项 目的 价值 重量 比率 关系 Pren 7 (Presis Piesizs Poesj s Pies) 表示 粒子 群 某 次 迭代 的 最 优 解 , 全 局 
只 有 四 种 情况 ， 利 用 这 一 特性 提出 了 GR-DKP(Greedy Repair 极 值 Ses = (es Biz Be t Brea) 表示 种 群 在 进化 过 程 中 的 全 局 最 
Algorithm for D{0-1}KP) 算 法 ， 设 计 了 基于 GR-DKP 的 粒子 群 优 解 ， 个 体 极 值 和 全 局 极 值 由 适应 度 函 数 确定 ， 粒 子 进 化 公式 为 
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一 求解 算法 (the PSO based Greedy Repair Algorithm for solving vj mV +C R (Pher X) Ca T (gen X) (4) 
D(0-1)KP, PSO-GRDKP)U!, 文献 [3] 在 研究 D(0-1) KP. 第 一 遗 1 
传 算法 (First Genetic Algorithm，FirEGA) 中 再 次 提出 贪心 修复 sende lee 6) 
与 优化 算法 (Greedy Repair and Optimization Algorithm, "m n if n > sig(vi?) (6) 
GROA)PBI, GROA 算法 与 GR-DKP 算法 均 按 照 非 递增 项 价值 d otherwise 


密度 对 项 进行 贪心 选取 , 若 项 集中 有 多 项 是 选取 状态 时 ,选择 其 中 : tR vi RRE t RAIRE A E 
价值 密度 最 大 项 。 杨 洋 等 人 进一步 优化 了 GROA 算法 ， 构 造 子 的 第 j 维 速度 ， 沐 表示 在 1 次 进化 迭代 中 第 i 个 粒子 的 第 j 
了 新 贪心 修复 优化 算法 (New Greedy Repair and Optimization 维 位 置 。 (4) 式 为 粒子 群 的 进化 动力 方程 , 主要 方面 构成 : 

Algorithm，NGROA)I， 该 算法 也 按照 非 递 增 项 价值 密度 对 项 节 属 于 粒子 个 体 的 惯性 势 ，c (Cov 一 六 ) 来 自 粒 子 个 体 当前 
进行 贪心 选取 , 但 当 项 集中 有 多 项 是 选取 状态 时 , 选择 项 集 的 历史 最 好 位 置 D 的 引力 势 ， 代 表 粒 子 “ 个 体 认 知 ”; 
介 值 最 大 项 ， 应 用 NGROA 算法 提出 NFirEGA(New First c hes - 蕊 ) 是 来 自 群 体 当前 历史 最 好 位 置 ees 的 引力 势 , 代 
Genetic Algorithm) 的 研究 结果 表明 较 之 FirEGA 提升 了 D1{0- 表 粒 子 “ 社 会 认 知 ”09。 学 习 因 子 c,c 作为 粒子 “个 体 认 知 ” 
KP 求解 质量 。 最近， 文献 [8] 定 义 了 项 集 价值 密度 概念 ， 提 和 “社会 认 知 ”的 主要 加 权 系 数 ， 表 示 进 化 过 程 对 粒子 个 体 
按 非 递增 项 集 价值 密度 对 项 进行 贪心 选取 的 策略 , 当 项 集中 信息 和 社会 信息 的 学 习 继承 程度 ， 主 要 影响 着 粒子 的 优化 
项 是 选取 状态 时 ， 选 择 满足 解约 束 条 件 的 价值 密度 最 大 项 。 标识 别 能 力 。 7i,5,5 表 示 (0, 1) 之 间 的 随机 数 。 
述 贪 栖 修复 与 优化 算法 B53 的 时 间 复 杂 度 为 000 。 多 数 研究 PSO-GRDKP 算法 外 是 D{0-1}KP 的 一 种 粒子 群 求解 算法 ， 
采用 文献 [3] 提 出 的 GROA 算法 设计 D{0-1}KP 的 各 类 启发 ”算法 搜索 空间 维度 D=3n-1， 粒子 适应 度 函 数 对 应 个 体 解 的 背 

式 进化 算法 ， 如 差分 进化 算法 (differential evolution algorithm， 包装 入 价值 ,学 习 因 子 ae 的 取 值 是 2, 算法 时 间 复 杂 度 是 O07), 
DE)9 、 变 异 蝙 晤 算法 (Mutated Double Codes Binary Bat 该 算法 利用 贪 禁 修 复 与 优化 算 子 对 种 群 中 不 可 行 粒子 进行 校正 

Algorithm ，MDBBA)U0、 差 分 进化 帝王 蝶 优 化 启发 式 算 法 和 优化 ， 有 效 提升 种 群 中 优质 个 体 比 率 ， 增 强 算法 寻 优 能 力 。 
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(Monarch Butterfly Optimization with Differential Evolution , 1.2 不 可 行 个 体 贪 梦 修复 与 优化 算法 
DEMBO)0I、 飞 蛾 搜索 算法 (Moth Search Algorithm, MS)?! 、 何 妆 朝 等 人 在 研究 D{0-1}KP 问题 时 先后 提出 GROADI] 、 


基于 Lagrange 插值 的 学 习 猴 群 算法 (Lagrange Interpolation — GR-DKPIJ 两 种 贪 禁 修复 算法 , 其 中 GR-DKP 在 研究 粒子 群 


based Learning Monkey Algorithm，LSTMA)U3、 基 于 环 论 的 演 算法 (PSO-GRDKP) 时 提出 ， 以 下 给 出 GR-DKP 算法 伪 代 码 。 
化 算法 (Ring Theory Based Evolutionary Algorithm，RTEAI)U4、 算法 1  GR-DKP 
基于 离散 混合 教学 的 优化 算法 (Discrete Hybrid Teaching 输入 : 3n 个 项 的 价值 向 量 P 、 重 量 向 量 w ; 背包 容量 C So Hons: 
Learning based Optimization Algorithm，HTLBO)035 等 ， 以 上 研 按 非 递增 项 价值 密度 次 序 保存 各 项 下 标 ; Aoo: 待 修复 粒子 。 
究 常 选用 FirEGA 算法 中 和 基本 启发 式 算 法 作为 基线 ， 通 过 对 输出 ;修复 优化 后 的 可 行 粒子 mo. 和 适应 度 值 fo 。 
比 求解 质量 和 收敛 速度 以 论证 所 研究 算法 的 优化 性 能 。 1. fweight =0, fralue=0 
在 采用 第 一 数学 模型 和 GROA 算法 构造 D{0-1}KP 启发 2. FOR i IN range(0.n-l) ; 
REAR ARPA, SCRRDIO]IBS D FESCUSESIRISH], TE IF O; Higi mua 2D: 
UDKP 实例 和 SDKP 实例 上 , 双重 编码 二 进 制 蝙蝠 算法 (Double Xai = 0, X3; = 0, X32 =0 
codes Binary Bat Algorithm, DBBA) 的 求解 质量 比 FirEGA 差 ， 7 
在 IDKP 实 例 和 WDKP 实例 上 ,DBBA 的 求解 质量 好 于 FirEGA:; xj = fei = wj 
帝王 蝶 优 化 算法 (Monarch Butter fly Optimization, MBO) 在 best, 3.  k=3n-1 
mean 及 worst 三 个 评价 指标 上 的 求解 结果 要 明显 差 于 FirEGA 4. WHILE weism>c AND £20: 
JEU! PSO-GRDKP 算法 求解 精度 和 稳定 性 优 于 FirEGATUI, IF( m=! ): 
这 些 研究 结果 表明 ,采用 同一 类 贪 禁 修复 与 优化 方法 的 不 同 启 fveghi- = wgi] xpi] = 0 
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5. FOR J IN nee(o3n-D : 
e(l) 
IF( soma mia = 0 AND feit wan < © ): 
X, = 1, fweight* = w, 
6. FOR i IN rnge(0,3n-D): fvalue+ — x, x p, 
7. RETURN x,fvalue 
算法 1 58 2-4 步 是 对 粒子 x 的 修复 操作 、 第 5 步 是 对 粒子 x 的 优 
化 操作 , 算法 时 间 复 杂 度 是 0(n) 。 将 算法 1 第 2 P emasa?) 
置换 为 j=arg max usassn (Pi) ， 则 算法 演化 为 NGROAD。 


2 ”基于 项 集 的 D{0-1}KP 粒子 群 优化 算 ; 


2.1 基于 项 集 的 粒子 修复 与 优化 算法 贪 禁 策 略 
D{0-1}KP 已 有 研究 提出 的 不 可 行 粒 子 贪 焚 修 复 策 略 
B41， 均 按照 项 价值 密度 大 小 对 数据 排序 预 处 理 。 本 文 引入 


项 集 价值 密度 R(ie{0,1…,n-1D ,按照 R 大 小 以 项 集 为 单位 对 
数据 进行 非 递 增 排列 ， 之 后 按 次 序 选取 项 集 的 项 。 文 献 [8] 给 
出 了 三 种 玉 算 法 如 下 : 

R = max(e,€,,4. €, ) (7) 
R? =X esu (8) 

S pi 
R= ta (9) 

È wiu 


易 见 ，R! 算法 等 效 于 GR-DKPIM 策 略 ， 结 合 NGROA 算 
AE, EX RTF: 


RÈ =e (10) 

以 下 给 出 项 集 ; 两 种 项 选择 策略 item 下 : 
item} = max (e; |x; =1, j € {3i,3i 13i  2)] (11) 

— 3i 2, pir x =1 
Hemp Le lx, -Lje(3i3ix1)), dtf (12) 


合式 (7)~(10) 和 式 (11)(12), 得 到 八 种 D{0-1}KP 基于 项 集 
的 贪心 修复 优化 方法 (Group Greedy Repair and Optimization 
Algorithm, GGROA)， 见 表 1 所 示 。 
表 1 D{0-1}KP 粒子 群 算法 的 GGROA 策略 


Tab. 1 GGROA of D {0-1} KP for PSO algorithm 
项 集 价值 密度 Rie{0.…,n-DD R m RO R 
项 集 i 的 项 选择 策略 item; item? 


GGROA 的 伪 代 码 描述 见 算法 2， 参数 me{12.3, 各 对 应 式 
(7)~(10) 的 四 种 项 集 价值 密度 ， 参 数 ze{1, 弛 对 应 式 (11)(12) 的 
两 种 项 选择 策略 。 

算法 2 GGROA( Hr",item: ) 

输入 : 3n 个 项 的 价值 向 量 P 、 重 量 向 量 w ;背包 容量 Cc ; Bu: 
按 非 递增 项 价值 密度 次 序 保存 各 项 下 标 ; tossed: 待 修复 粒子 。 
输出 :修复 优化 后 的 可 行 粒子 son 和 适应 度 值 CD) 。 

1: FOR ;IN »ange(0,n—1l) ; 


Flag; 20 
2* fweight = 0, fvalue = 0 
3: FOR i IN range0.n-D) ; 


IF Qux X3 2D: 
23i —0,X3;44 = 0, X3j42 =0 
j = item? 
xj = 1, fweight+ = wj, Flag; =1 
4: k=n—l 
5: WHILE eign>c AND >0: 
IF( Fagy =! ): 


| 
j= item; 
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weight = wjXy — 0, Flags -0 
大 = 大 一 1 


6: | FOR j IN wse0.n-D): 


，_ HUE 
了 = item; 


IF( Flag ymp] — 0 AND weight +w; < C ) : 


x, =l, fweight+ = w, 


N 


FOR i IN wgse(0,3n-D): fvalue+ = x, x p, 
8: RETURN x, fvalue 
算法 2 第 2 步 ~ 第 5 是 对 粒子 x 的 修复 操作 、 第 6 步 是 
对 粒子 x 的 优化 操作 ， 算 法 时 间 复 杂 度 是 O00 。 
2.2 基于 GGROA 的 D(0-1)KP 粒子 群 优化 算法 
文献 [和 在 构造 PSO-GRDKP 算法 中 , 采用 了 基础 离散 粒 
子 进 化 式 (4)~(6)， 学 习 因子 的 取 值 ee 是 2。 为 避免 粒子 进化 
公式 和 进化 参数 对 研究 内 容 的 影响 , 方便 同一 基准 条 件 下 同类 
研究 间 的 辨析 ， 以 下 采用 文献 [4 和] 的 方案 构造 PSO- 
GGRDKP(PSO based GGROA for solving D{0-1}KP)， 伪 代码 
见 算法 3， 其 中 参数 mz 算法 2。 
算法 3 PSO-GGRDKP(?".:) 
输入 : o 项 的 价值 向 量 P 、 重 量 向 量 w ,背包 容量 C ,种 群 规模 w ， 
进化 代数 了 ， 学 习 常 数 =e =2 
输出 :最 佳 粒子 sw 和 装 入 背包 总 价值 的 近似 最 优 解 finesse o 
1: FOR ;IN range(0,n-l) ; 
R” = function, (P..W,) 
排列 R ERmetoes-0, POWS 
E FEN Hin 
3: 初始 化 种 群 : population - (x, v) |i ea.) 
4: FOR i IN "enge. N) : 


价值 密度 R 非 递增 次 序 将 项 


poi 


uM 


x,, fitness (x, ) -GGROA(H" ,item: ) 
5 best = index _ max ($ fitness (x, ) | i e {1,--, N3]) 
6 Ebest = Phest = Xpest 
7: r=0 
8 WHILE ż <T: 
9 FOR i IN ranged.N) ; 
10: FOR Jj IN nnse(03n-D : 

v? EV EC n TX) +C, n, x) 


IF (x > sigQ7): x? =0 


ELSE: x'-! 
11: xil, fitness (x ) -GGROA(H" i tem ) 
12: IF Cfimess (x) > fitness (pps )) : p, =x" 
13: IF (fimess (Pres ) > fitness (8 pest )): 8 = Pru 
14: t=t+1 


15: RETURN gyes, fitness (8 pest) 
算法 3 (P nnn Æ (0, 1) 区 间 的 随机 数 ， 第 1 步 中 
Junction, (m € (12.3,4]) 对 应 四 种 项 集 价值 密度 计算 函数 ， 即 式 
(7)~10)， 第 10 步 sig(x) 采用 式 (5)， 算 法 时 间 复 杂 度 是 
O(nlogn)+20(n)+20(nN)+Tx[O(nN)+0(n)] , HF N«nHT«n, 
itt PSO-GGRDKP 算法 时 间 复 杂 度 是 OG) 。 


3 ”仿真 实验 与 结果 分 析 


3.1 实验 方案 设计 

D1{0-1}KP 数据 集 是 观察 和 评测 算法 性 能 的 标准 数据 集 B]， 
数据 集 由 逆 强 相关 D{0-1}KP 实例 IDKP)、 强 相关 D(0-1) KP 
实例 (SDKP)、 弱 相关 D{0-1}KP 实例 (WDKP) 和 不 相关 D (0- 
1}KP 实例 (UDKP) 四 类 数据 组 成 ， 四 类 D{0-1}KP 实例 的 数 
据 规 模 分 别 为 300<3n<3000 。 为 验证 D{0-1}KP 的 各 类 贪 禁 修 
复 策略 性 能 特点 ， 进 行 两 组 实验 。 实 验 一 : 采用 八 种 基于 项 
集 贪 焚 修 复 优 化 策略 的 D{0-1}KP 粒子 群 算法 实验 ， 以 探究 
八 种 基于 项 集 贪 禁 修复 优化 策略 的 可 行 性 、 性 能 以 及 与 数据 
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实例 类 型 的 适应 性 关系 。 实 验 二 : 评测 PSO-GRDKPÜ!, PSO- 表 2 D {0-1}KP 的 GGROA 策略 及 序号 
NGROADKP PSO-GGRDKP 三 种 典型 不 可 行 个 体 贪 禁 修 Tab.2  Ggroas and their serial numbers of D{0-1} KP 
复方 法 的 D(0-1) KP 粒子 群 性 能 特点 。 两 组 实验 均 以 比较 id (R item) id (Ritem) 
求解 时 间 和 解 计 算 结果 分 析 算 法 性 能 。 所 有 实验 均 在 1 (1,1) 5 Q.D) 
ThinkStation P330 计算 机 上 进行 , 电脑 配置 Intel® Core™ 2 (1,2) 6 2) 
17-8700 CPU-3.2GHZ., 16 GB DDR4 Jj f£. NVIDIA P2200 3 (2,1) 7 (4,1) 
显卡 ， 操 作 系 统 是 Microsoft Windows 10 教育 版 ， 算 法 4 (22) 8 (42) 


均 使 用 Python3.6 编码 。 表 3 所 示 为 数据 规模 3n=900 和 3n=1800 的 四 类 实例 上 运 
id D{0-1}KP 实例 通过 基本 动态 规划 法 计算 得 出 的 最 优 。 行 八 种 贪 栖 修复 优化 方法 的 粒子 群 算法 实验 数据 ，PSO- 
TEC opt 、 粒 子 群 算法 的 粒子 规模 为 200、 进 化 代数 同 项 集 妆 GGRDKP 算法 名 称 下 标 编 号 对 应 表 2。 


tT 


n ， 粒 子 群 算法 独立 运算 20 次 的 近似 最 优 解 最 大 值 为 best. 应 用 八 种 贪 禁 修复 策略 对 八 个 数据 实例 IDKP3、 

平均 值 为 mean、 最 差 值 为 worst、 平 均 时 长 为 7 。 IDKP6、SDKP3、SDKP6、UDKP3、UDKP6、WDKP3、 

3.2 实验 一 算法 数据 与 分 析 WDKP6 分 别 独立 求解 20 次 ,20 个 求解 结果 分 布 情况 箱 
表 2 列 出 八 种 项 集 贪心 修复 策略 的 对 应 编号 (id)。 线 图 如 图 1 所 示 。 


表 3 D1{0-1}KP 标准 数据 实例 PSO-GGRDKP 算法 实验 结果 
Tab.3 PSO-GGRDKP algorithm experimental results of D (0-1) KP standard dataset 


实例 名 PSO- PSO- PSO- PSO- PSO- PSO- PSO- PSO- 
(最 优 值 ) GGRDKP: GGRDKP; GGRDKP; GGRDKP4 GGRDKP; GGRDKPs GGRDKP; GGRDKPs 
best 234772 234737 233983 234160 234450 234413 234720 234668 
iens mean 234565.65 234592.45 233706 233832.9 234296.05 234139.9 234490.8 234354.6 
worst 233790 234183 233092 232980 234103 233119 234204 232898 
pug T(s) 12.18 9.07 11.03 11.72 10.88 10.06 11.22 9.40 
ERR 0.001 0.001 0.005 0.004 0.002 0.003 0.001 0.002 
best 452250 451982 450537 450171 450971 450816 451829 451876 
IDKP6 mean 451303.05 450804 449725.95 449256.7 450148.5 449898.55 451071.3 450849.7 
(452463) worst 449531 448999 446422 447609 448009 447226 449593 446822 
T(s) 46.70 36.34 48.26 41.80 46.64 42.04 44.23 3747 
ERR 0.006 0.008 0.013 0.011 0.010 0.012 0.006 0.012 
best 237313 237638 237273 237428 237414 236682 236973 236671 
SDKP3 mean 236644.85 236362.7 236428.45 236276.9 236661.5 235866.1 236064.6 236051.5 
(238248) worst 235198 232226 233959 232525 234957 234252 234342 234145 
T(s) 11.48 9.56 10.46 9.10 10.66 9.01 10.37 9.32 
ERR 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 
best 462883 461926 462498 460133 462797 462300 462197 462010 
ESI mean 461454.95 458730.3 460320.15 456405.1 461389.9 459623.25 460858.45 460292.2 
(466097) worst 459209 451076 456090 447800 458299 454888 455773 457203 
T(s) 45.41 37.77 41.34 36.29 41.70 35.94 41.77 37.06 
ERR 0.01 0.02 0.01 0.02 0.01 0.01 0.01 0.01 
best 256286 255916 255712 255529 255879 255844 256083 256061 
mean 255516.15 255491.25 255003.15 254928.75 255034 255448.45 255300.9 255713.75 
dn worst 253410 254768 252796 253758 253558 254668 253481 254466 
T(s) 10.77 8.98 12.47 9.00 10.50 10.32 11.77 9.85 
ERR 0.004 0.004 0.006 0.007 0.006 0.005 0.005 0.004 
best 464720 464210 463503 463134 463414 463434 464072 464274 
mean 463532.4 462831.4 462020.7 461409.2 461794.25 462586.1 462883.55 463214.55 
ADS worst 462308 459786 459081 459144 459866 460557 461176 461365 
(465030) T(s) 44.32 38.93 44.48 3644 45.07 40.76 43.48 37.56 
ERR 0.005 0.007 0.009 0.010 0.009 0.007 0.007 0.006 
best 263944 267162 262269 267172 265006 267221 264710 267102 
UDKP3 mean 256631.2 266185 252478.55 266165.65 259174 266483.15 258556 266184.5 
(184006) worst 227502 263431 229837 263871 233367 264297 237289 263701 
T(s) 11.83 10.50 12.01 10.57 12.11 10.59 11.77 10.66 
ERR 0.05 0.01 0.06 0.01 0.04 0.01 0.04 0.01 
best 520183 529834 517776 530221 517264 529030 521267 528480 
UDKP6 mean 506575.15 526417.6 504295.95 526795.8 504110.5 526440.7 508569.7 526602.7 
(536578) worst 457476 519485 464104 522972 451696 521096 464333 524425 
T(s) 47.47 41.07 47.79 41.75 47.41 41.24 47.59 41.31 
ERR 0.06 0.02 0.06 0.02 0.06 0.02 0.05 0.02 


aats 


由 表 3 和 图 1 可 以 看 出 ， 粒 子 群 算法 采用 八 种 项 集 贪 禁 ”结合 每 种 数据 实例 的 平均 时 长 7 ,进一步 选择 各 类 数据 实例 时 
修复 策略 都 是 可 行 的 , 但 不 同 的 项 集 贪 焚 修 复 策 略 在 同一 类 数 ， 间 性 能 较 优 的 求解 算法 , 其 中 ,IDKP: PSO-GGRDKP;,, SDKP: 
据 实例 上 存在 显著 性 能 差异 , 相同 的 项 集 贪 焚 修 复 策略 在 不 同 ^ PSO-GGRDKPs, WDKP: PSO-GGRDKPs, UDKP: PSO- 
类 型 实例 上 性 能 表现 也 有 不 同 。 同类 数据 的 两 个 实例 , 确定 习 GGRDKP2. 

均 最 优 近似 解 前 三 的 贪心 修复 策略 交集 为 该 类 数据 较 佳 的 算 3.3 实验 二 算法 数据 与 分 析 

法 ， 其 中 ，IDKP: PSO-GGRDKP;, PSO-GGRDKP», SDKP: PSO-GRDKP、PSO-NGROADKP、PSO-GGRDKP 算法 实验 
PSO-GGRDKP;, PSO-GGRDKPs, WDKP: PSO-GGRDKP:. 数据 见 表 4 所 示 , 其 中 PSO-GGRDKP 算法 根据 实验 一 分 析 结 
PSO-GGRDKPs, UDKP: PSO-GGRDKP», PSO-GGRDKPe. ”选择 四 类 数据 实例 上 的 性 能 较 优 的 贪 禁 修复 策略 进行 实验 。 
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-r:IOBHgIT.| TT I IR. 
(e) WDKP3 (£WDKP6 (g)UDKP3 (h)UDKP6 
图 1 D {0-1}KP 实例 PSO-GGRDKP 求解 结果 箱 线 图 
Fig. 1 Box plots ofPSO-GGRDKP results for D [0-1] KP instances 
表 4  PSO-GGRDKP, PSO-NGROADKP, PSO-GRDKP 算法 实验 结果 
a Tab.4 Algorithm experimental results of PSO-GGRDKP, PSO-NGROADKP, PSO-GRDKP 
PSO-GGRDKP PSO-NGROADKP PSO-GRDKP 
Dataset C opt 
best mean worst T/s best mean worst T/s best mean worst T/s 
IDKP1 61500 70106 70090 70057 70037 1.106 70098 70097 70090 1.430 70106 70105 70098 1.284 
IDKP2 103936 118268 118232 118168 118037 4.272 118235 118235 118232 5.294 118268 118268 118268 4.862 
IDKP3 214453 234804 234640 234326 232619 9.216 234785 234778 234679 11.666 234802 234795 234686 11.048 
IDKP4 251980 282591 282565 282350 281660 16.444 282579 282570 282423 20.327 282571 282562 282442 19.352 
IDKP5 297482 335584 335460 335077 334317 29.298 335580 335562 335383 31.347 335541 335538 335474 30.543 
IDKP6 415217 . 452463 . 451933 450909 449351 47.024 452410 452359 451760 47.535 452458 452412 451920 47.303 
IDKP7 434677 489149 488657 487982 486639 50.667 489044 488975 488469 65.012 489118 489095 488870 58.752 
IDKP8 464860 . 533841 532893 531437 527685 84.127 533749 533604 533060 88.726 533833 533810 533398 75.925 
IDKP9 454989 528144 527458 526463 524077 92.545 528010 527978 527678 104.692 528115 528106 527965 95.705 
IDKP10 496541 581244 580528 578794 576617 111.433 581043 580968 580350 142.916 581207 581194 581020 118.073 
SDKP1 60143 94459 94325 93875 93098 1.178 94431 94411 94258 1.792 94232 93962 93413 1.423 
SDKP2 100200 160805 160211 159881 159009 4.652 160663 160631 160359 5.340 160671 160321 159616 5442 
SDKP3 151680 238248 236996 236248 234059 12.609 238132 238081 237774 13.362 238055 237476 236394 11.994 
SDKP4 214644 340027 338126 337472 336426 18.292 339817 339562 338477 21.282 339786 339055 337662 20.905 
SDKP5 300632 . 463033 459165 457548 453168 28.383 462667 462163 461013 32.062 461309 459479 456502 33.365 
SDKP6 289412 . 466097 . 462003 460618 457818 41.004 465596 465323 464135 49.359 465194 464176 462190 47.159 
SDKP7 396421 620446 — 616378 612841 605344 56.047 0619691 618860 617739 65.980 0618204 616231 612674 64.845 
SDKP8 424240 670697 . 664952 662384 654486 73.405 0669310 668958 667365 83.399 0668735 667212 664795 83.565 
SDKP9 465564 739121 734123 730826 724267 93.362 737018 736819 735649 105.648 736476 734747 731842 107.098 
SDKP10 477837 765317 . 758980 753068 747259 124.712 762795 762668 761977 126.717 761900 760266 757622 132.407 
WDKP1 65423 83098 83025 82932 82731 1.122 83085 83082 83051 1.336 | 83083 83067  Á 82927 1.361 
WDKP2 102795 138215 137843 137759 137494 4.354 138175 138163 138015 5.128 138213 138185 137928 5.356 
WDKP3 196603 256616 . 256114 255722 255289 9.431 256470 256432 256054 11.409 256586 256369 255310 12.101 
WDKP4 239770 315657 314849 314640 314193 16.692 315585 315521 315022 20.221 315620 315418 314548 21.317 
WDKPS5 325793 428490 427124 42639] 424448 26.110 4428188 428068 427534 31.547 428424 428129 426928 34.034 
WDKP6 351265 466050 . 464687 463570 461558 37.787 465756 465690 465292 46.473 465857 465599 464261 50.446 
WDKP7 411219 547683 545903 544410 542740 51.114 547239 547156 546451 63.677 547369 547091 545921 70.765 
WDKP8 424552 576959 574803 572924 569842 67.464 576408 576351 575937 86.628 576677 576236 574893 92.975 
WDKP9 478991 650660 647825 645650 643211 85.902 650031 649728 648581 112.974 650271 649897 648347 110.999 
WDKP10 493578 678967 676082 674452 670718 109.833 678113 677977 677273 137.185 678742 678413 677273 140.024 
UDKP1 56429 85740 85690 85576 84893 1.282 85686 85667 85354 1.371 85676 84357 78486 1.383 
UDKP2 94596 163744 163173 162263 159583 4.829 163645 163513 162327 7.232 163744 162729 156376 5.388 
UDKP3 184006 269393 267595 266626 263915 10.756 268885 268721 267718 11.982 269200 266037 252792 12.004 
UDKP4 204182 347599 344978 343223 337971 18.971 346991 346867 345717 25.439 347336 344506 333362 21.206 
UDKP5 279856 442644 . 435563 432351 426821 29.374 440623 440151 438481 34.637 440938 434662 415259 33.438 
UDKP6 344578 536578 529459 525276 522693 42.012 534097 532875 531044 47.067 534091 526749 503799 48.258 
UDKP7 424973 0635860 631620 626296 616479 57.390 0632897 631970 629597 61.907 632642 625851 600343 66.428 
UDKP8 412859 650206 . 637976 634314 626440 76.533 646661 646379 644234 84.020 645782 639408 618604 87.228 
UDKP9 437234 718532 705861 700529 694344 96.546 713920 713702 712034 104.043 712556 705909 684026 113.5 
UDKPIO 463681 779460 | 762985 754952 745142 120.374 772963 772437 770927 133.303 771741 766303 746945 130.33 
从 表 4 看 出 ， 三 类 贪 禁 修 复 优化 策略 的 粒子 群 算法 贫 禁 修复 优化 粒子 群 算法 四 类 数据 实例 的 平均 解 误差 


求解 质量 和 时 间 性 能 各 不 相同 。 定 义 粒 子 群 算法 的 解 误 
差 率 SRR=1-”%m ， 各 数据 实例 的 解 误 差 率 平均 值 为 
ERR_AVE ， 该 值 越 小 ， 算 法 性 能 越 优 。 进 一 步 统计 三 类 


率 与 算法 平均 时 长 见 表 5 所 示 。 表 中 三 类 贪 禁 修 复 优 化 
粒子 群 算法 分 别 是 PSO-GGRDKP, PSO-GRDKP, PSO- 
NGROADKP 。 


a 
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d 5 看 出 ， 在 四 类 数据 实例 上 ， 三 类 粒子 群 算法 的 3)PSO-GGRDKP, PSO-NGROADKP, PSO-GRDKP 算法 
均 解 误差 率 呈 现 一 致 趋势 : IDKP<WDKP<SDKP<UDKP , 即 对 解 误 差 率 与 数据 实例 相关 系数 、 项 价值 密度 均值 、 项 价值 密 
IDKP 实例 求解 误差 率 最 低 , 而 UDKP 实例 的 解 误差 率 最 高 。 度 标准 差 等 数据 特征 高 度 相 关 。 

三 类 算法 中 ，PSO-NGROADKP 平均 解 误差 率 最 低 ， 仅 为 以 上 结果 表明 ，PSO-GGRDKP 算法 的 解 误 差 率 略 高 于 
0.2% ， PSO-GGRDKP 的 平均 解 误差 率 略 高 于 PSO- PSO-NGROADKP、PSO-GRDKP, fH £i 5 E USE T D(0-1) KP 


iN 


T 


le 


NGROADKP, PSO-GRDKP 分 别 是 0.7%、0.4%。 在 四 类 数据 的 算法 时 间 性 能 。 

实例 上 ，PSO-GGRDKP 的 时 间 性 能 均 显著 优 于 PSO- 本 文 以 粒子 群 算法 为 例 ， 构 造 D(0-1)KP 的 PSO- 
GRDKP 和 PSO-NGROADKP，PSO-GGRDKP 算法 时 间 GGRDKP 算法 并 验证 其 性 能 。 受 各 类 应 用 GROA 的 改进 启 
性 能 较 之 PSO-GRDKP 提升 12.9% ， 较 之 PSO- 发 式 算法 研究 结果 启示 ， 本 文 推断 在 不 同 启发 式 求解 算法 中 
NGROADKP 算法 时 间 性 能 提升 13.8%。 三 类 粒子 群 算 使 用 GGROA 的 求解 性 能 可 能 会 有 差异 ， 限 于 篇 幅 ， 这 将 作 
法 的 四 类 数据 实例 时 间 性 能 趋势 表现 不 一 致 ，PSO- 为 进一步 研究 工作 讨论 。 近 期 ， 文献 [20] 提 出 D{0-1}KP 问题 


GGRDKP、PSO-NGROADKP、PSO-GRDKP St 1X tT [Ra] 44 粒子 群 算法 的 项 集 个 体 编 码 方案 ， 接 下 来 考虑 结合 文献 [20] 
能 表现 最 好 的 数据 实例 各 自 为 :WDKP、SDKP 和 IDKP。 的 工作 改进 PSO-GGRDKP 算法 。 
RS 三 类 粒子 群 算法 平均 解 误差 率 与 平均 时 长 


Ht 


Tab. 5 Solution average error rate and average time of three kinds of PSO 参考 文献 i 
" . . PSO-GGRDKP . PSO-NGROADKP PSO-GRDKP [1] Guder J. Discounted knapsack problems for pairs of items [D]. 
TUR ERR AVE T AVE ERR AVE T AVE ERR AVE T AVE Nuremberg: University of Erlangen-Nurnberg, 2005. 
IDKP 0.002 44.613 0.000 51.895 0.000 46.285 D] Guldan B. Heuristic and exact algorithms for discounted knapsack 
SDKP 0.010 45.364 0.002 50.494 0.005 50.817 problems [D]. Nuremberg: University of Erlangen-Nuremberg, 2007: 1- 
WDKP 0.005 40.981 0.001 51.658 0.001 53.938 78. 
UDKP 0.017 45.807 0.005 51.100 0.015 51.916 B] A, EER, FLR, X. 基于 遗传 算法 求解 折扣 {0-1} 背 包 
AVE 0.009 — 44.191 0.002 51.287 0.005 . 50.739 题 的 研究 [J]. ea 2016, 39 (12): 2614-2630. (He Yichao, 
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性 ， 表 6 给 出 了 四 类 数据 实例 的 项 价值 系数 与 项 重量 系 discounted (0- 1! knapsack problem [J]. China Journal of Computer, 2016, 
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